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The Schrodinger equation is solved for the case of a particle confined to a small region of a 
box with infinite walls. If walls of the well are moved, then, due to an effective quantum nonlocal 
interaction with the boundary, even though the particle is nowhere near the walls, it will be affected. 
It is shown that this force apart from a minus sign is equal to the expectation value of the gradient 
of the quantum potential for vanishing at the walls boundary condition. Variation of this force with 
. . . time is studied. A selection of Bohmian trajectories of the confined particle is also computed. 
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^ ; I. INTRODUCTION 

' f-H '■ In most of the problems of quantum mechanics the Hamiltonian of the system is time dependent and so one needs to 
O I solve the time dependent Schrodinger equation (TDSE). Problems with moving boundary conditions are an interesting 
^ I ^ , class of such time dependent problem. Such a system was first considered by Fermi [l| in connection with the study 
■ of cosmic radiation. After, several authors studied problems with moving boundaries [il-Q i @ ^^'i references therein. 
[ Different aspects of the problem of a particle in a one-dimensional infinite square-well potential with one wall in 
^ ■ uniform motion have been discussed by earlier authors. Exact solution of the time dependent Schrodinger equation 
I ^ i ' for this problem at first, as far as we know, was given by Doescher and Rice [2[. Schlitt and Stutz Q considered the 
application of the sudden approximation to the rapid expansion of the well. Pinder investigated the applicability 
CO ■ of both adiabatic and sudden approximation for both expanding and contracting wells: it was shown that sudden 
^ , approximation is appropriate to the expanding well provided the wall speed is sufficiently great, but this approxi- 
. mation may not be applied for the contracting well irrespective of the rate of contraction. Using the semiclassical 
approximation, Luz and Cheng Q evaluated the exact propagator of the problem. The energy gain and the transition 
amplitudes and probabilities between initial and final energy eigenstates of the problem have been calculated 0] . A 
recent numerical study of a particle in a box with different laws for the movement of the wall show that physical 
quantities like probability density and expectation value of position or mean value of the energy have a smooth be- 
havior for small speed of the moving wall. In contrast, if this speed becomes large, many irregularities appears as 
sharp bumps on the probability distribution or a chaotic shape on the averaged values of position and energy . 

The aim of the present paper is to probe some aspects of the time dependent boundary condition for a particle 
confined in an infinite square well that have remained hitherto unnoticed. 

Let us focus on the effect of the time dependent boundary condition while one of the infinite boundary walls in a 
^ ' box is moved where we have a well-localized Gaussian wave packet which remains peaked at the center of the box, 
5^ ■ Xc, well away from walls. Now, by calculating the effective quantum force [l^ one can study the way this effective 
quantum force changes with time. Due to such a force, the expectation value of the momentum in the direction 
perpendicular to the walls gradually changes in time. Then one can compare curves of the quantum effective force 
for the static (when the wall is at rest) and the dynamic (when the wall moves) situations, and pinpoint the instant 
from which the dynamic curve deviates from the static one. Such an instant shows the time at which the confined 
particle begins to feel the motion of the wall. 

Physically one expects that when the width of the initial Gaussian packet is much smaller than the initial width of 
the box, (To ^ ^Oi one might consider a Gaussian wavepacket to be realizable in the box trap. But, this approximation 
casts some doubts upon the computations (i.e. is the effect an artifact of the tails of the Gaussian at the boundaries?). 
So for this purpose, in addition, we consider a localized state with a finite support embedded in the support of the 
expanding box trap at time t = 0. It is most natural to consider the particle-in-a-box eigenstates of a tiny box 
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centered at x = Xc, and suddenly released at t = to become the initial state without the necessary approximations 
to assume a Gaussian packet as an initial state. Such state is called "tiny-box state" afterwards. 

The computed Bohmian trajectories [lll - [T3| for the static and the dynamic situations are also instructive in revealing 
the conceptual ramifications of such an example. In Bohm's model each individual particle is assumed to have a 
definite position, irrespective of any measurement. The pre- measured value of position is revealed when an actual 
measurement is done. Over an ensemble of particles having the same wave function ip, these ontological positions are 
distributed according to the probability density p = j-i/jp where the wave function ijj evolves with time according to the 
Schrodinger equation and the equation of motion of any individual particle is determined by the guidance equation 
V = j / p, where v is the Bohmian velocity of the particle and j is the probability current density. Solving the guidance 
equation one gets the trajectory of the particle. 

The plan of this paper is as follows. Section |TT] contains a very brief review of the relevant mathematical steps 
leading to the exact solution for the problem. In Section IIIII numerical computations related to the effect of the time 
dependent boundary condition are presented. Finally, in Section [IVI we present the concluding remarks. 



II. BASIC EQUATIONS 

Consider a narrow box inside a wide box with a particle inside the inner one. Walls of the outer box are at = 
and xji = £o and waUs of the inner one are at xi = (^o — ^i)/2 and X2 = {£o +^i)/2 where £i <^ £o- At time t = the 
inner box is suddenly removed and the right wall of the outer box starts to move uniformly with velocity u. Infinite 
wall speed u corresponds to a hard wall at a; = 0. We discuss solution of TDSE for two cases. Initial wavefunction to be: 



i) a Gaussian wave packet well localised in the center of the tiny box. In fact in this case we have a truncated 
Gaussian packet, because of confinement of the Gaussian packet with infinite tails in a narrow region, and so 
the name "truncated Gaussian packet" (TGP)). The problem concerning the tails of Gaussian packet that was 
mentioned in the introduction, now is translated to the truncation. 



ii) the ground state of the narrow box with kick momentum k (tiny-box state, "TBS" for abbreviation). 
To get a picture see fig. (1). 

Using the propagator of a rigid box with the left wall at x^ = and the moving right wall in a constant velocity u 
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(1) 



and the relation, 

= J dx'K{x,t; x\0)iJoix') , (2) 
one gets the wavefunction at any time having ^'0(2;) in hand. £{t) = £0 + ut shows the position of the moving wall at 



time t. At this stage, it must be mentioned that due to the Galilean invariance of the Schrodinger equation |13j . the 
case of both moving walls is equivalent to the case of one wall in motion but with u as the relative velocity of walls. 
With the initial wavefunction to be a TGP well localised in the center of the box Xc = ^o/2. 
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where Q{x) is the step function; one gets. 
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where 
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and Erf is the error function: Erf(z) = e * dt. In the second case we take the initial wave function as, 

Mx) = y^sin[^(a;-xi)] e'^(^-"=) e(a;-xi) e(.T2-x) . 



(5) 



(6) 



In this case the relation of il:{x,t) is cumbersome; there are eight modified error functions, Erfi(z) = e* dt in 

its summand. 

In original Bohm approach to causal interpretation of quantum mechanics [ill [isj to introduce the concept of 
particle, Schrodinger equation is decomposed into two real equations by expressing the wavefunction in polar form 
■0 = Re^^^^. Then, vector filed v = p/m is constructed from the vector filed p = VS* and assuming that v defines at 
each space-time point the tangent to a possible particle trajectory passing through that point. In this interpretation 
of quantum mechanics one gets, 



dp 

'dt 



-V{V + Q) 



(7) 



where Q = — (7i^/2m)V^i?/i? is known as quantum potential. Analogous to classical physics, in Bohm's model of 
quantum theory one has, 



dt 



'dt' 



(8) 



where the mean value is defined for an ensemble of density E? and momentum p = VS*. In the standard approach 
to quantum mechanics the right hand side of eq. (jS]) is meaningless (ref. [l^ pp: 111-113). Using eq. ^ and taking 
the expectation value of eq. ([7]) one obtains, 



dt 



-(v(y + g)). 



(9) 



For the case of a particle within a box with one wall moving, using the integration by part one can find 
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where we have used the fact that wavefunction is zero on both walls. The general case of boundary condition will be 
considered in the appendix. One obtains, 
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where the second equaUty holds at the boundaries only. In consequence we have, 



d(p) _ 
dt ^ ^2m 

= .fqm{t) J 
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(10) 



where /qin(i) was called quantum effective force by Dodonov and Andreata [To]. It must be mentioned that in the 
context of standard approach to quantum mechanics one can obtain eq. IjlOp by simultaneous application of the 
Schrodinger equation, 



2m dx^ 

and time-derivative of the expectation value of momentum operator 

d{p) 
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as it was done at first by Dodonov and Andreata for the case of an impenetrable wall at a; = 0. 
Now, taking the integral of both sides of cq. ((TO)) leads to 
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where in the case of TGP (p)(0) = hk Erf[^^^] whereas for the case of TBS (p)(0) = hk. 



III. NUMERICAL CALCULATIONS 



In this section we work in a unit system where 7i = 1 and m = 0.5. Other parameters are chosen as £i ~ £o/20, 

Co = ^i/lO and £o = 1. Conservation of the probability, given by J^^^^ dx = 1, can be used as a parameter 
that gives us a test on the precision of the results (It must be noted in the case of TGP because of truncation total 
probability is not equal to unity, instead it is equal to Erf[ ^J^^ ] — 0.9999994267). We use Simpson's rule for taking 
the integral of eq. Using the conservation of the probability as control variable, we got a good numerical stability 
by just taking the first 400 terms of infinite sum appearing in the relation of ip^x^t) in both TGP and TBS cases. 
Using the RungeKutta method for solving the guidance differential equation a selection of Bohmian trajectories is 
presented. 

Fig. [T]show the initial wave function for both TGP and TBS and the walls of the well to give a feeling about how 
narrowly the initial wavefunction is. 

In fig. [3] a selection of Bohm paths is presented for a initially motionless wavefunction, i.e., fc = in cqs. ([3]) and 
From parts a) and c) one can see that Bohm particle in the static case remains at rest at the center of the box. 
We have checked long time behaviour of Bohm trajectories that starts at the tail of the leading half in the dynamic 
situation and saw that for TGP, particle eventually moves on a path approximately parallel (approximately, because 
it has very small oscillation around the parallel path) to the path of the moving wall, i.e., with the velocity of wall, 
but for TBS it moves with a velocity less than the velocity of moving wall. 

Fig. [3] show the expectation value of position operator for a initially motionless wavefunction. From this figure 
and fig. [2] one finds that in the static case, i.e., fixed wall, the Bohm path which initially placed at the centre of the 
motionless packet, xq ~ {x){0) (in the case of TBS xq ~ 0.5 whereas xq ~ (l/2)Erf[5/\/2] = 0.4999997 in the case 
of truncated Gaussian), moves with the centre point subsequently, x{t) = {x){t), at least up to time 0.003 that we 
have considered. But, this is not true in the dynamic case. The reason is the quantum effective force which has been 
shown in fig. |4l In static case fqm{t) is zero for both tiny-box and truncated Gaussian states. Deviation of {x){t) 
from its static value (x)(0) take place sooner for the TBS. This shows that the particle begins to feel the motion of 
the wall sooner for the case of TBS compared to the case of TGP. In figures [5] and El we have plotted time-dependence 
of quantum effective force for a fixed value of the speed of the moving wall, u = IOOtt, but different values of kick 
momentum k. Noting these figures one finds that for a moving packet, i.e., fc ^ 0, quantum effective force is not 
zero in the static case contrary to the case of a motionless one. Comparison of these figures show that particle begins 
to feel the motion of the wall approximately twice sooner in the case of TBS in comparison to truncated Gaussian 
packet. 
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FIG. 1: (Color online) Initial probability density. Vertical blue dashed lines show the walls of the narrow box and the vertical 
solid blue lines stand for the walls of the wide box. 



Fig. [6] reveals that in the case of truncated Gaussian packet, quantum effective force is the same for both static and 
dynamic cases for A: < and it is zero for fc > in the dynamic case, at least for our parameters and time-domain 
t G [0,0.0005]. /qm(0 deviates from zero in positive direction for fc < but in negative direction for fc > 0, i.e., 
particle accelerates for fc < but decelerates for fc > 0. 

Noting fig. [7]which displays /qm(0 for the motionless TBS and TGP but for different wall's speed, one finds: 1) in 
the presented region of time, quantum effective force for TGP is negligible compared to the TBS (one must note that 
in longer time limit opposite behaviour take places according to fig. SJ, 2) direction of deviation from zero changes 
with u and 3) deviation time increases with u. 



IV. SUMMARY AND DISCUSSION 



In this paper we studied the solution of TDSE for a particle in a) tiny-box state and b) truncated Gaussian 
packet of an infinite square well with one wall in uniform motion. We showed that due to a quantum effective force, 
which apart from a minus sign is the expectation value of the gradient of the quantum potential in the context of 
Bohmian mechanics, the expectation value of the momentum operator changes gradually with time. We studied the 
variation of this quantum effective force with time for different values of the speed of the moving wall in the case of 
a motionless packet and different values of the kick momentum but fix value of the speed of the moving wall. Some 
Bohm trajectories for the motionless packet were also plotted. We have learned from the numerical calculations that 
the particle in TBS begins to feel the motion of the wall sooner in comparison to TGP. This may be understood by 
computing the speed of propagation [l^ for both TBS and TGP. Other ramifications of this study, like a contracting 
box, dependence of quantum effective force on related parameters like mass of the confined particle, width of the 
initial packet, other initial packets like excited particle-in-a-box eigenstates and other types of boundary conditions, 
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FIG. 2: (Color online) A selection of Bohm trajectories, x{t), for it — lOOvr and an initially motionless wavefunction: a) 
TBS, static case; b) TBS, dynamic case; c) TGP, static case and d) TGP, dynamic case. In each figure black curve starts at 
xo — Xc — 2(70, red one at xq = Xc and the green one at xo = Xc + 2ao- 



like periodic ones, call for further consideration. 
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FIG. 3: (Color online) Expectation value of position operator versus time, {x){t), for u — IOOtt and an initially motionless 
wavefunction: a) TBS and b) TGP. In each figure the black curve shows {x}(t) for the dynamic case, the green curve shows 
{x){t) for the static case and the red one shows the Bohm trajectory for the dynamic case which starts at 2:0 = (a;)(0). 



Appendix A: General form of quantum effective force for a particle in a box 



1. Bohmian Mechanics 



In the context of Bohmian mechanics, the actual momentum of the particle is given hy p = dS/dx, where S/h is 
the phase of the wave function. Thus, for the time-derivative of the expectation value of the actual momentum of the 
confined particle inside the one-dimensional box, one has [l^ 
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where we have used the Leibniz's formula, and lit) — d£{t)/dt. Now, using the continuity equation 
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and the generalized Hamilton- Jacobi equation 
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FIG. 4: (Color online) Quantum effective force versus time, /qm(i), for an initially motionless wavefunction: a) TBS and b) 
TGP. In each figure the black curve shows /qm(<) for the dynamic case and the red one shows /qm(i) for the static case. 



which can be obtained by putting the polar form ^{x,t) ~ R{x,t)e^^^^'*'>/^ in the Schrodinger equation, one gets 



Integrating by part, leads to 
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Thus, using the fact that inside the box the classical potential is zero, one has 
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Now, using the definition of quantum potential Q = —{h'^/2m)V'^R/R, 
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FIG. 5: (Color online) Quantum effective force versus time, /qm(i), for the tiny-box state and u = lOOvr: a) k = —75n, b) 
k = —50tv, c) k = —25n, d) k — 25tv, e) k = 50n, f) fc = 75n. In each figure the black curve shows /qm(t) for the static case 
and the red one shows /qm(t) for the dynamic case. 



Last integral can be evaluated as follows, 



dx = 2 



dx^ R ) dx 
Finally, one obtains 
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for the quantum effective force. Now, one can use the general form (jA6|) for other types of boundary conditions like 
periodic ones, ^/'(0,t) = ■ip{£{t),t) and il;'{Q,t) = ip'{£{t),t), for which the momentum operator does have eigen-statcs 
obeying this periodic boundary condition. 



2. Standard Quantum Mechanics 

Simultaneous application of equations (|11[) and (jl2p with Leibniz's formula leads to 
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FIG. 6: (Color online) Quantum effective force versus time, /qm(t), for the truncated Gaussian packet and u = IOOtt: a) 
k — ~75tv, h) k = — SOtt, c) fc = —25n, d) k = 25n, e) k = 50n, f) fc = 75n. In each figure the dotted black curve shows /qm(t) 
for the static case and the red one shows /qm(t) for the dynamic case. 



(p) is real, so, one can write 

where Re(2;) and Ini(2) shows the real and imaginary part of z, respectively. Eq. (|A8I ) is converted to eq. (|A6|) by 
using the polar form of the wave function. It is easily seen for "vanishing at the walls" boundary condition eq. (|A8[) 
yields the simple form (|10p . 




dtp 
dx 
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FIG. 7: (Color online) Quantum eflFective force versus time, /qm(t), for a motionless a) TBS and b) TGP in dynamic case (force 
is zero in static case for fc = 0). In each figure black curve is for u — 207r, red one is for it = lOOvr and the green one is for 
u = 200n. 



